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We consider a family of spherical three dimensional spacelike slices embedded in the Schwarzschild 
solution. The mean curvature is constant on each slice but can change from slice to slice. We give 
a simple expression for an everywhere positive lapse and thus we show how to construct foliations. 
There is a barrier preventing the mean curvature from becoming large, and we show how to avoid 
this so as to construct a foliation where the mean curvature runs all the way from zero to infinity. 
No foliation exists where the mean curvature goes from minus to plus infinity. There are slicings, 
however, where each slice passes through the bifurcation sphere R — 2M and the lapse only vanishes 
at this one point, and is positive everywhere else, while the mean curvature does run from minus to 
plus infinity. Symmetric foliations of the extended Schwarzschild spacetime degenerate at a critical 
point, where we show that the lapse function exponentially approaches zero. 
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I. INTRODUCTION 

Constant mean curvature (CMC) foliations of the 
Schwarzschild geometry have been constructed by Brill, 
Cavalho and Isenberg These foliations degenerate 
when the lapse "collapses"; foliations in the vicinity of 
this critical point have been been investigated in In 
both of these papers, it is assumed that the trace of the 
extrinsic curvature is not only a constant on each slice, 
but retains this constant value from slice to slice. In this 
paper we analyse more general families of CMC folia- 
tions, with the trace K of the extrinsic curvature that 
can change with time. The "collapse" of the lapse can 
again be described by analytic approximations. 

The standard way of viewing General Relativiy as a 
dynamical system is by considering the 4-manifold as fo- 
liated by a sequence of spacelike 3-surfaces [!, Q • Each 
3-surface inherits a 3-metric, gij, and an extrinsic curva- 
ture Kij = (l/2)C n gij [j| where £ n is the Lie derivative 
along the normal, and thus K^ is a geometric object, the 
analogue of the time derivative of the metric. Each slic- 
ing is equivalent to a choice of time function. A standard 
choice is to demand that the trace of the extrinsic cur- 
vature, g lJ Kij, usually written as K and known as the 
mean curvature of the surface, be constant along each 
slice. Hence these are the constant mean curvature, or 
CMC, slices. CMC slices are attractive for a number of 
reasons. The Einstein equations give an evolution equa- 
tion for K, it is (in vacuum) 

BK 

V 2 N~K^K lJ N = NC n K=—-N l d l K, (1) 



and this can be regarded as an equation for the lapse 
function, N, of a CMC slicing. This is a nice linear ellip- 
tic equation that satisfies the maximum principle. 

In the standard, conformal, methods of constructing 
initial data for the gravitational field, choosing the trace 
of the extrinsic curvature to be a constant simplifies the 
equations. One gets a single, nonlinear scalar equation 
for the conformal factor instead of a coupled system of 
nonlinear equations [6|]. 

By taking the trace of the equation defining Kij , one 
can show C n y/g = ^fgK. This tells us that K is just 
the fractional time rate of change of the volume along 
the normal, and, in a cosmology, the CMC slices are the 
'Hubble time' slices, with the instantaneous value of K 
equalling the Hubble 'constant'. In Minkowski space, on 
the other hand, the 'mass hyperboloids', t 2 — r 2 = m 2 , 
are CMC slices, with K = 3/m. Slices in general asymp- 
totically flat spacetimes mimic the behaviour of the mass 
hyperboloids, in that they are everywhere spacelike but 
become null at null infinity. This feature makes them of 
value for the analysis of gravitational or any other form 
of radiation, see e.g. Q [!]• 

One major advantage of CMC slicings for numerical 
relativity is that, if we consider gravitational waves of 
a fixed wavelength, the number of wave cycles that a 
CMC slice intercepts is finite. To resolve a wave in a 
numerical computation, we need the separation between 
data points to be less than the wavelength. This means 
that the domain of a code can extend all the way to null 
infinity, and track the waves all the way out with a finite 
number of grid points. This phenomenon is independent 
of whether one compactifies or not. 



where N and N z are the lapse and shift, respectively. 

If K is a spatial constant, the shift term drops out and 
the equation reduces to 



f)K 

V 2 N - K v Ki . N= ™, 



(2) 



II. EXPLICIT CMC FOLIATIONS 

We possess a great deal of understanding about the 
spherical CMC slices of the Schwarzschild solution (2j. 
We can write the 3-metric and extrinsic curvature ana- 
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lytically in terms of the Schwarzschild radius R. They 
depend only on two parameters, K, the trace of the ex- 
trinsic curvature, and C, that quantifies the transverse- 
traceless part of the extrinsic curvature. The expressions 
are 
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The extrinsic curvature is diagonal, with the non-zero 
components being 
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plays a key role. It has several meanings, k = dR/dL, 
where L is the proper distance along the slice, therefore 
2k/R\s the 2- mean- curvature of the round 2-spheres as 
embedded in the 3-slice. Also k — Njc, where Nk is the 
Killing lapse, the dot product of the 'timelike' Killing 
vector with the unit normal to the slice. For small values 
of K and C, k 2 is positive for both small and large R, 
with only two roots. We must have positive k 2 for the 
metric to make physical sense. Let us label the larger of 
these two roots Rt- One can show easily that R t < 2M. 
The zone where k 2 is positive from R = R t to R = oo 
corresponds to a slice that starts out at one future null 
infinity (if K > 0), reaches a minimum surface at R = 
Rt (hence l t' for 'throat') and continues on to the other 
future null infinity. If C < 8M 3 K/3 it crosses below the 
bifurcation point, if C > 8M 3 K/3 it crosses through the 
upper quadrant. As C becomes larger, while holding K 
fixed, the polynomial k 2 rises up, the two roots approach 
each other until they finally touch at a value of R which 
we call R*. The value of K and C at that point equal 



2R+ - 3M 
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This is a critical point of the slicing, As (C, K) change 
so as to approach this point one gets the 'collapse of the 
lapse' and 'slice stretching' that one is used to in maximal 
slicing^, K we increase C even further, k 2 is posi- 
tive over the entire range R = (0, oo). When this is true, 
the CMC slice starts at null infinity and plunges into the 
singularity. We have previously analysed the CMC slices 
where we hold K fixed and just allowed C to change Q. 
Obviously there exists a much richer class of CMC slic- 
ings, where both C and K change. One really only needs 
to give some relationship between C and K. However it 
is easier to introduce some parameter time t and write 
C(t) and Kit). One can always reparametrise, which 
will change the form of both C(t) and K(t). However, 
the ratio, (dK j dt) / (dC / dt) , is unchanged. 
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This equation is solved by the lapse function N, where 



N = (3k + k 



dr- 
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Here /3 is a (time-dependent) constant. One can verify 
that the lapse equation ((2|), is solved by N as given in 
Eq.©. 

It is clear that Eq.@ is a linear elliptic inhomoge- 
neous equation, so therefore the solution can be written 
as a linear combination of a 'particular' solution of the 
inhomogeneous equation combined with a solution of the 
homogeneous equation (V 2 — K l i Kij)tp = 0. Since k is 
the Killing lapse, we know that 



(V 2 - K ij K i:j )k = 0. 



(9) 



We also know that we can change C without changing K. 
This also gives us a lapse that satisfies the homogeneous 
equation, call it <p(C). It can be chosen to depend linearly 
on C (one may have to subtract off some multiple of 
k). Finally, the inhomogeneous solution can be can be 
adjusted so as to linearly depend on K (again, one may 
have to subtract off multiples of k and C). This is why we 
get the three constants /?, C, and K appearing explicitly 
in Eq.©. 

Let us now restrict ourselves to the case where C and K 
are small enough that k has a zero at R = R t . Let us also 
assume K > 0. This means that the slice comes in from 
one future null infinity, has a minimal surface at R = Rt 
and goes out to the other future null infinity. We know 
that the Killing lapse k, the solution to Eq.©, is anti- 
symmetric, passing through zero at R = Rt- By adding 
an appropriate multiple of k we can find the solution that 
is proportional to K and the solution proportional to C 
that are symmetric around the throat. 

Let us have a spherical CMC slicing of the 
Schwarzschild solution. This can be described as a curve 
in (C, K) space. There is still some freedom. The term 
in the lapse proportional to k can be thought of as 'pure 
gauge'. It slides the slice along the Killing vector with- 
out changing either K or C. Since we are interested in 
the effects of changing K and C we would like to elimi- 
nate this freedom. The obvious way is to adjust the lapse 
by adding an appropriate amount of k so that the lapse 
is symmetric around the throat. This gives us a Neu- 
mann boundary condition dN/dL = dR/dLdN/dR = 
kdN/dR — at R = R t where L is the proper distance 
along the slice. Since k = at the throat, this looks like 
a trivial equation. It is not, because dN/dR generically 
blows up at the throat like l/k. This has the effect of 
reducing the number of constants in the solution from 
three to two. The second freedom is the choice of time 
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parametrisation of the slicing. It appears that the 'natu- 
ral' choice is to set (3 = 1 in Eq.©. This allows the time 
translation vector of the slicing to coincide with the time- 
like Killing vector at null infinity. It turns out that we 
can satisfy both conditions, i.e., (3 = 1 and the symmetry 
at the throat, simultaneously. 

It is interesting to look at expression Eq. © in the limit 
of large R. In that limit k » KR/3. Therefore the 
integrand is dominated by —9K/K 3 r 2 . This integrates 
to 9K/K 3 r, so the integral equals -9K / K 3 R+0(1/ R 2 ). 
The k outside the integral becomes KR/3 so we finally 
get 



3K 

N~(3k- — + 0(1/R) 



(10) 



Therefore N/k — > (3 for large R. Since k is the Killing 
lapse, if P = 1 we can have the time translation of the 
slicing equal the Killing vector. With this choice the time 
parameter is the retarded time. 

While expression Eq. © is easy to understand near in- 
finity, it is difficult to see from it what happens near the 
throat. To find the cleanest expression, we have to ma- 
nipulate Eq.©. It is easy to see that 
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and the expression for the lapse can be rewritten as 
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Now integrate by parts to get 
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If we change the range of integration from the throat to 
R, instead of from R to oo we get 
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It turns out that if we choose 
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i.e., eliminate the term that is of the form constant xk we 
get the symmetric lapse function that satisfies dN/dL = 
at the throat. This is 
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We obviously assume that we are in the subcritical regime 
so that R t exists. It is straightforward, if tedious, to show 
that this expression satisfies the CMC equation, Eq.([2]). 

As stated earlier, we can simultaneously set dN/dL = 
at the throat (just by choosing the lapse as given by 
Eq. (fT7|) ) and simultaneously set (3 = 1. We start off with 
a curve in (C, K) space on which we put some coordinate 
time label. This defines C and K . Now compute 



(18) 
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and scale the time by 7, i.e., 

dt = "/dt. 



(19) 



This rescales C and K by 7. Now the value of (3 for which 
the two lapse functions, Eq.© and Eq.fTT]), are equal is 
= 1! 

The quantity M + K 2 r 3 /9 + KC/3 - 2C 2 /r 3 that ap- 
pears twice in the expression for N is nothing but (r 2 /2 
times) the first derivative of k 2 and this is positive on 
the entire interval [Rt, 00). It only goes to zero as we ap- 
proach the critical point. It is the 'going to zero' of this 
function that gives the classical 'collapse of the lapse'. 

This lapse, while maintaining the CMCness of the 
slices, will not keep the metric in the form Eq. ©, 
because the normal to the slice is not along the R — 
constant direction. It is easy to work out what the ap- 
propriate shift is since the Killing vector is along this 
direction. We know that the Killing lapse is olk = k and 
we also know that a 2 K — (3 2 K = 1 — 2M/R. Therefore 
(3 2 K = {KR/3 - C/R 2 ) 2 => (3% = k(-KR/3 + C/R 2 ). In 
turn we get N R = N{-KR/3+C/R 2 ). Since g RR = 1/k 2 
we get N R = N(-KR/3 + C/R 2 ). This argument works 
for both choice of lapse, Eq.© and Eq. p"?"]) . In the up- 
per half plane, if N > 0, we want N R < for large R 
because the normal to the CMC slice is leaning over more 
towards null infinity than the Killing vector is. This shift 
does not vanish at the throat for the symmetric N as dis- 
tinct from the Killing shift. It will be positive (if N is 
positive) there and we seek a symmetric solution. This 
means that we have a discontinuity in the shift across 
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the throat. This is to be expected, because, as we move 
forward in time, the coordinate range of R expands. The 
opposite will occur in the lower half plane. From the for- 
mula for the symmetric N, Eq. (fTT)) . we can easily read 
off the value of the lapse at the throat as 



N{R = R t ) = 
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This is obviously positive if dC/dt - dK/dt(Rf/3) > 0. 
If we have a symmetric slicing and if N — at any one 
point then it must be zero on an entire closed 2-sphere. If 
dK/dt is positive, we can immediately see from applying 
the maximum principle to the CMC equation, Eq.([2]), 
that if the lapse is positive in the center it is positive 
everywhere. One way to show this is to assume the op- 
posite. Let us assume that we have a region on the 3- 
surface where N > in the interior, and N = on the 
boundary. Therefore dN/dR < on the boundary. Mul- 
tiply the lapse equation, Eq.([2]), by N and integrate over 
the region in which TV > 0. One gets an immediate con- 
tradiction! Therefore we have a foliation instead of just 
a slicing. This kind of argument was first introduced by 
fiol ] . This expression for the central lapse in Eq. (fP7j) gives 
us an immediate consistency check. Let us consider the 
situation where the variations in C and K are such that 
the radius of the throat does not change. Since R = Rt is 
the zero of k 2 , by inspecting the formula for k 2 , Eq.©, 
the condition is that 
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But, of course, if the radius does not change the lapse at 
the throat must vanish. This is exactly what we get from 
Eq.®- 



and K to change with time, i.e., either Eq.® or Eq. (fT7]) . 
This has the form 
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IV. MATCHING CMC TO MAXIMAL SLICING 

Using the explicit form of the metric(s) in the previ- 
ous section, we can show that one can smoothly match a 
spherical CMC slicing (with either a constant or time- 
dependent K) to a maximal slicing in an extended 
Schwarzschild solution. The matching is performed along 
a timelike surface with fixed Schwarzschild radius R — 
Rj > 2M. We can handle the case where the CMC slices 
are on the 'inside' and the maximal slices are 'outside', 
and the converse. 

We wish to match a patch of spacetime with metric 



with 



/ N 2 (a 2M\ N I C 

W \W 3 I ' w 

o o 



¥) 






R 2 











R 2 sin 2 9 J 



N 
~k 



dr- 



C 



3 



r 2 fc 3 



(22) 



III. TWO SPACETIME METRICS 

Let us start with a CMC slice whose 3-metric is given 
byO and let us drag this along the timelike Killing vector, 
i.e., we set K = 0, C = 0. We can now write down the 
4-metric associated with this slicing (it is not a foliation 
because the lapse vanishes at the throat). It is 
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along some cylinder of constant Schwarzschild radius R 
Rj to a patch of spacetime with metric 
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Another spacetime metric can be constructed using the 
same 3-metric but choosing a lapse function that allows C 



W i th fc2 = 1 _2M + ^ 

We need to check that the Israel-Darmois junction 
conditions @ are satisfied. Note that we allow (3 to 
be two different time dependent functions, one in each 
patch. Really, we only need /3_fc(t) to change with time, 
it makes sense to keep 0^ = 1. This guarantees that the 
time in the 'maximal' zone is the proper time at infinity. 
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The Israel-Darmois junction conditions are that intrinsic 
metric and the extrinsic curvature of the 3-dimensional 
matching surface be continuous. We can actually arrange 
that the whole 4-metric be continuous. 

The key condition is that we choose C, K and C to 

satisfy C — C — —j 1 - ■ This has to be satisfied on every 

~ • KR 3 

time slice. Therefore we also want C — C This 

guarantees that k = k along the matching surface. We 
also can use the free parameter /3k (t) to maintain N = 
N. This guarantees that the 4- metrics to the left and 
right of the surface R = Rj are the same. 

Now we need to look at the extrinsic curvature. The 



condition C — C 



KR 



- guarantees that Kgg = Kgg and 



K<f„f, — K (/,(/,. This means that we only need to check Kqq- 
The normals to the surface R — Rj are respectively — 
(0, i.0, and = (0, p0, 0. Hence = n M . We have 

^oo = «o;o = -^ob n R = -^g Rli (2gon,o- 9oo,n)n R . Since 
the metrics match on the surface, we have <7o/j,o = 9o^,o- 
Therefore the problem reduces to comparing {N/k) r to 
(N/k), R . Looking at Eqns.© and ||23J). This reduces 
~ • KR 3 

to comparing C to C ^ J -. These are obviously equal. 

Therefore we can match both the intrinsic metric and the 
extrinsic curvature along the 3-surface defined by R = 
Rj. 

The metric along the matching surface is C°° while the 
matching of the extrinsic curvatures guarantees that the 
metric perpendicular to the surface is C . Further, con- 
sider the Hamiltonian and momentum constraints along 
the matching surface. The behaviour of the metric and 
extrinsic curvature means that there is no jump discon- 
tinuity in ( 3 ) R ~ K ab K ab - K 2 or V a (K ab - g ab K) across 
the surface. This means that there are no stresses along 
the matching surface. 

Consider the Hamiltonian and momentum constraints 
on a spacelike CMC slice. These are (3) i? - K ab K ab + 
K 2 = and \7 a (K ab - g ab K) = 0. If we write 
K ab = K? + l/3Kg ab we can rewrite the constraints 



as WR-KfKl b 



fore the initial data can be regarded as maximal data 
coupled to a constant negative-density 'matter' field at 
rest with 1 67173 = -(2/3)K 2 . Therefore we can regard 
this as a negative cosmological constant. We can even 
naturally deal with the situation where the the cosmo- 
logical constant is time dependent since we are free to 
allow K to change from slice to slice. 

Thus we can consider this matching as gluing a sphere 
in a maximally sliced anti De Sitter spacetime to a max- 
imal Schwarzschild exterior. This would involve holding 
K fixed in time. Alternatively, we could have the inte- 
rior maximal and the exterior be CMC. Obviously, we are 
only considering the Schwarzschild spacetime in a spher- 
ically symmetric slicing. This has no gravitational radi- 
ation. Nevertheless, the community of people who nu- 
merically analyse binary black hole collisions might find 
such slicings appealing. The black holes could be put in 
the maximal zone while the radiation could be analysed 



in the asymptotically null CMC zone. 



V. FOLIATIONS VERSUS SLICINGS 

If dK/dt is positive and if dC/dt - dK/dt(R%/3) > 0, 
then the lapse function is positive everywhere and we 
have a foliation. Why is this? How restrictive a condition 
is it? 

If we hold K fixed and increase C we automatically 
get a foliation Q. The lapse is everywhere positive and 
the slices move forward in time. However, if we hold C 
fixed and increase K, the slices move forward near both 
infinities and backwards in the middle. One way of seeing 
this is to look at the (implicit) formula for R t , 
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Now compute dRt/dK. We get 
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This shows that R t increases as K increases. However, 
in the upper quadrant of the Schwarzschild solution the 
Schwarzschild radius decreases as one moves forward in 
time from the bifurcation sphere (at R = 2M) to the 
singularity at R = 0. 
Similarly, we have 
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Therefor R t decreases with increasing C and so the slice 
moves forward in time. 

Therefore to have the slices moving forward every- 
where, when K changes, one needs to simultaneously 
increase C as well as K. This behaviour is illustrated 
in the two figures. 

One way of finding CMC slices in Schwarzschild is 
by means of a height function approach, see 0, [Io| . 
This involves considering a slice defined as t — where 
t = T — h(R), and where T and R are the standard 
Schwarzschild coordinates. The condition that the slice 
be CMC reduces to a second order differential equation 
for h. This can be explicitly integrated once to give 0] 
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FIG. 1: K changes but C is fixed, slices can intersect. 




integral, just as we do! 

We can get an interesting slicing that is 'almost' a foli- 
ation by setting C = 8M 3 K/3. All these slices have their 
throats at the bifurcation sphere, R = 2M. These slices 
all touch at R = 2M, and so the lapse is zero there. Oth- 
erwise, as K increases, the lapse is positive everywhere 
else. This slicing allows K to run all the way from — oo 
to +00. If we use K as our label time, the lapse function 
of this slicing is 

1 8M 3 - i? 3 
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This slicing covers all of the left and right quadrants, but 
never enters the upper or lower quadrants, the 'black' 
or 'white' hole zones. This shows us that, while we can 
get a foliation that covers the range in K of [0, 00), we 
cannot find a foliation that allows K to run the whole 
range (—00, +00) because the Schwarzschild solution is 
time symmetric and a foliation that gets to K = +00 
must break this. To some degree, this is a word game. 
If you prespecify the range to be covered, i.e., one seeks 
a foliation that goes from K = -D/M to K = +D/M, 
where D is a large number, one can do this. Start off with 
a moment of time symmetry slice , i.e., (K = 0, C = 0) 
and choose a curve in (K, C) space so that dC jdK is 
only infinitesmally bigger than i? 3 /3. This will reach any 
desired value of if, in particular K = D/M. Now add the 
time reversal of this, and we have the desired object. Of 
course, all such foliations eventually run into the critical 
curve at some finite value of K. This, however, would be 
bigger than the specified D/M. 



VI. CRITICAL FOLIATIONS 



FIG. 2: The trace K is constant and C changes; the slices 
form a foliation. 



In retrospect, the fact that the lapse function can be 
written in such a compact form should not have been 
too surprising. The lapse function can be viewed as the 
derivative of the height function with respect to K and 
C. Therefore one might expect to write the first derivative 
of the lapse as a function, and the lapse itself as a simple 



Henceforth we assume the subcritical regime — a min- 
imal surface at Rt and the Neumann boundary condition 
dN/dL = where L is the proper distance along the slice. 
This guarantees that the slices are symmetric about the 
throat. The line element, expressed in terms of comoving 
time t and the areal radius R, takes the following form 

2m" 
~R~ 



ds 1 



-dt 2 1 



A - 



c_ 

2A«^ 



KR 



-dtdR + 



die 

k 2 



+ R 2 dfl 2 



(30) 



Here A 
grals 

X 



1 + dr C r -±li K ■ Let us define following inte- 



lK 2 r 2 



R fc (2M+^ + 22p 



IC' 2 



rdr 



Y = 



r 2 ( 



6m + 2KC 



2AC' 



(2M- 



2Kr 



2KC 
3 



4C 2 



7 dr. (31) 
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The formula for the lapse function, N 
rewritten as 



fcA, can be 



N 



R±dK_ , dC 
3 dt ~*~ dt 



M 



K 2 R 3 , jfC 
9 ' 3 

2dK,„ n dC lxr 
--—kY + 2—kX. 
3 dt dt 



2C 2 
R3 



(32) 



Define T 



\B\k 



1/2 



The asymptotic behaviour of the 



lapse function near the critical point follows from (|33[) . 
the analysis of the decaying of k 2 near R* and ([57)1 . One 
obtains 



(38) 



This representation of the lapse N is convenient when 
there exist minimal surfaces. From (|32[) one clearly sees 
that if the minimal surface exists at Rt, k(R t ) = 0, then 



N = 2- 



-K + C 



M 



K 2 R? 



KC 
3 



2( 



(33) 



One finds from (j3"2")) that at the throat dN/dL = 
hm k2 ( : + 2< ^ X _ ■ Therefore the condition 

dN/dL = yields the differential equation 



C 



-1 

2X 



3X ' 



(34) 



this is a highly nonlinear relation, since both integrals X 
and Y depend in a convoluted way on C and K . Never- 
theless, it is possible to give a compact analytic descrip- 
tion of this foliation near the critical point , AT* , of 
the foliation. At this point vanish both k 2 and the first 
derivative The integrals X and Y become divergent 



at the critical point, but the right hand side of (|34|) is 
finite everywhere and vanishes at the critical point. De- 
fine 



A = 



-K R t =- 

K 2 R 2 



3_C 2 

R 3 t 
c 2 



Ri 



12 



(35) 



Now assume that K, C, Rt are close to critical values 
A'* , C* , R« . One can show after a lengthy analysis and a 
number of careful estimates, that 



X 
Y 



i — Rt 
A 2 
-Rt 



'2k 



A 2 



2Rf 
A 



7T _ \ 72C 2 
2 ~ ) RtAni 



.(36) 



In the particular case of the trace K being independent of 
time, the estimate (|38|) coincides with the result derived 
earlier in Q. 



VII. CONCLUSIONS. 

We consider a family of spherical three dimensional 
spacelike slices embedded in the Schwarzschild solution. 
The mean curvature is constant on each slice but can 
change from slice to slice. One describes how the slices 
are stacked by defining the lapse function, that quanti- 
fies distance along the normal as one goes from slice to 
slice. We write down a simple expression for the lapse 
of any such slicing. This allows us to glue a patch of a 
Schwarzschild spacetime with a CMC slicing to a patch 
that is maximally sliced. It is easy to identify those slic- 
ings where the lapse is everywhere positive. The slices 
do not cross so one has a foliation. There is a barrier 
that prevents the mean curvature from becoming large, 
and we show how to avoid this so as to construct a fo- 
liation where the mean curvature runs all the way from 
zero to infinity. No foliation exists where the mean cur- 
vature goes from minus to plus infinity. However, if we 
consider the slicing where each slice passes through the 
bifurcation sphere, the point where R = 2M, we almost 
get a foliation because the lapse only vanishes at this one 
point, and is positive everywhere else, while the mean 
curvature does run from minus to plus infinity. There 
exist symmetric foliations of the extended Schwarzschild 
spacetime. They degenerate at a critical point. We show 
that the lapse function exponentially approaches zero at 
this critical point. 



The calculation is completely analogous to that of Sec- 
tion VII in and it is sketched in the Appendix. Define 
5= R t -R* and e = C* - C. Moreover, let SK = K t ~K* 
as R t tends to R*. Then (up to terms of 



and if 



lower order) A « 6k and e 



B X 2 



S 2 . Here A 



Hi 

2 



B = -2C, 
ical point 
equation 

Jt( e 



f- jAT^i?^ and is the value of k at the crit- 
The insertion of the above information into 
and (|3"6")) yields the differential equation 



1 



SKR: 



\B\^I 2 
V2i?? 



1 



-SKR: 



(37) 
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VIII. APPENDIX. 



The function k (R) vanishes at the minimal surface, at 
R t . Thus 2m + \KC = R t + \K 2 R^ + ^ and the de- 
nominator of the third factor of (|4l"j) can be represented 
as 



(A + iV) ■ (42) 
Here appears the function A, already defined in the main 
text. It is now convenient to replace r in the integral X 

byy = Jl-2±. Then X reads 



X = — 4y R t / dy 



F 



y/A + y 2 F 2 (A + y 2 F 3 



(43) 



Here F 4 = F z /y 2 for i = 2,3, and F 1 = Fx/ (l - y 2 ) . 
Now define a new variable z = y/yA; the integral (|43|) 
becomes 



Below we show how one calculates the integral X de- 
fined in the main text. Define following functions 



i 7 ! 

h 

Rt 



F, 







C 2 



K 2 r 



2„2 



r 

K 2 R 



2 D 2 



i r 2 
R? R 



9 y rij nt 
Rf\r r 2 

l K 2 R 3 t ( 4 - i 

9 4 \R$ 



R 3 t 



AC 2 



R t 



(39) 



The function k 2 defined in the main text can be written 
as follows 



k 2 



1 




(40) 



and the integral X takes the form 

dr 



X = -2JR t 



1 



l_«t VA^ 

r 



(2m + \KC + lK 2 r 3 - ^ 



(41) 



X 



A 2 



dzF\ (zy/A 



z 2 F 2 (zVE) (l + z 2 F 3 (zVA 

(44) 



One can split the integral J Q 
Jo 



1/VkS 



<-l/v'l0 4 Ki r-1 

Jo + h 



1/VkS 



into two parts: 
It is easy to check that 



the contribution coming from the second integral of the 
integrand of X goes to zero as 5 approaches zero. The 
first integral of the integrand of X in turn is well approx- 
imated by the integral 



X = -4- 



fRt 



-nR t 



dz 



1 + ^R t n{\ + z 2 R tK y 



(45) 



This can be explicitly calculated, with the result dis- 
played in the main text. Let us remark, that the function 
1/A explodes to infinity at the critical point of the folia- 
tion. Thus in this limit the quantity XA 2 becomes equal 
to — \72ft* and the first of equations ([36]) appears exact. 



The calculation of the other (Y) integral defined in Eq. 
H) proceeds in a similar way. 



